Power Law Scaling for a System of Interacting Units with Complex Internal Structure 
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We study the dynamics of a system composed of interacting units each with a complex internal 
structure comprising many subunits. We consider the case in which each subunit grows in a multi- 
plicative manner. We propose a model for such systems in which the interaction among the units is 
treated in a mean field approximation and the interaction among subunits is nonlinear. To test the 
model, we identify a large data base spanning 20 years, and find that the model correctly predicts 
a variety of empirical results. 
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In the physical sciences, power law scaling is usually 
associated with critical behavior, thus requiring a partic- 
ular set of parameter values. For example, in the Ising 
model there is a particular value of the strength of the 
interaction between the units composing the system that 
generates correlations extending throughout the entire 
system and leads to power law distributions [Q. In the 
social and biological sciences, there also appear exam- 
ples of power law distributions (incomes ||], city sizes 
extinction of species [Q, bird populations [||, heart 
dynamics |^). However, it is difficult to imagine that for 
all these diverse systems, the parameters controlling the 
dynamics spontaneously self-tune to their critical values. 

In this Letter, we raise an alternative mechanism by 
asking how power law distributions can emerge even in 
the absence of critical dynamics. The guiding principles 
for our approach, to be justified below, are: (i) the units 
composing the system have a complex evolving structure 
(e.g., the companies competing in an economy are com- 
posed of divisions, the cities in a country competing for 
the mobile population are composed of distinct neighbor- 
hoods, the population of some species living in a given 
ecosystem might be composed of groups living in different 
areas), and (ii) the size of the subunits composing each 
unit evolve according to a random multiplicative process. 

Fortunately, for one of the examples listed above, there 
is a wealth of quantitative data, and here we focus on 
a large database giving the time evolution of the size 
of companies Q]. In an economy, the units composing 
the entire system are the competing companies. In gen- 
eral, these companies have a complex internal structure, 
with each company composed of divisions (the subunits 
of each unit). It has been proposed that the evolution of 
a company's size is described by a random multiplicative 
process with variance independent of the size, and that 
each company can be viewed as a structureless unit 
However, later studies P-llq] reveal that the dynamics of 



real companies are not fully consistent with the simplified 
picture of Ref. Q . 

We develop here a model that dynamically builds a di- 
versified, multi-divisional structure, reproducing the fact 
that a typical company passes through a series of changes 
in organization, growing from a single-product, single- 
plant company, to a multi-divisional, multi-product com- 
pany . The model reproduces a number of empirical 
observations for a wide range of values of parameters and 
provides a possible explanation for the robustness of the 
empirical results. Due to our encouraging results for the 
case of company growth, our model may offer a generic 
approach to explain power law distributions in other com- 
plex systems. 

The model, illustrated in Fig. |l|, is defined as follows. 
A company is created with a single division, which has 
a size ^i{t — 0). The size of a company S = Yii£^i{t) at 
time t is the sum of the sizes of the divisions £,i{t) com- 
prising the company. We define a minimum size S'min 
below which a company would not be economically vi- 
able, due to the competition between companies; S'min is 
a characteristic of the industry in which the company op- 
erates. We assume that the size of each division i of the 
company evolves according to a random multiplicative 
process g. We define 
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where ?7i(t) is a Gaussian-distributed random variable 
with zero mean and standard deviation V independent 
of . The divisions evolve as follows: 

(i) If A^i(t) < S'min, division i evolves by changing its 
size, and £^i{t -|- 1) = ^^(t) + A£^i{t). If its size be- 
comes smaller than Smin — i-e. if ^i(t + 1) < Smin 
— then with probability pa , division i is "absorbed" 
by division 1. Thus, the parameter pa reflects the 
fact that when a division becomes very small it will 
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FIG. 1. Schematic representation of the time evolution of 
the size and structure of a company. We choose Smin ~ 2, and 
Pf = Pa = 1.0. The first column of full squares represents the 
size 5i of each division, and the second column represents the 
corresponding change in size A^i. Empty squares represent 
negative growth and full squares positive growth. We assume, 
for this example, that the company has initially one division 
of size ^1 = 25, represented by a 5 x 5 square. At t = 1, 
division 1 grows by A^i = 3. A new division, numbered 2, is 
created because A^i > Smin = 2, and the size of division 1 
remains unchanged, so for t = 2, the company has 2 divisions 
with sizes — 25 and ^2 = 3. Next, divisions and ^2 grow 
by 2 and —2, respectively. Division 2 is absorbed by division 
1, since otherwise its size would become ^2 = 3 — 2 = 1 which 
is smaller than 5111111. Thus, at time t = 3, the company has 
only one division with size ^1 = 25 + 2 + 1 = 28. Note that 
if division 1 would be absorbed, then division 2 would absorb 
division 1 and would then be renumbered 1. If, division 1 is 
absorbed and there are no more divisions left, the company 
"dies." 



to each company, and we find no significant dependence 
of the results on the value of £ ior £ — 20, 30 or 50. 
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FIG. 2. Probability density of the logarithm of company 
size for the model and for US publicly-traded manufactur- 
ing firms in the 1994 "Compustat" data base. To obtain 
these results, we assume that Smin is drawn from a trun- 
cated log-normal distribution with average value 5 x 10^. 
The numerical simulations were performed with parameters 
V = 0.15, Pf = 0.8, Pa = 0.05, and ^ = 50 (for these pa- 
rameter values, the actual probability of a new division being 
created per division and per iteration is approximately 0.01). 

It is common to study the logarithm of the one-year 
growth rate, ri = lni?i, where Ri = S{y + 1)/S{y), and 
S{y) and S{y + 1) are the sizes of the company in the 
years y and y + 1. The empirical distribution of ri for 
companies with size S is consistent with an exponential 
form llSfl 



no longer be viable due to the competition between 
companies. 

(ii) If A^j(i) > S'inin, then with probability (1 — p/), we 
set ^i{t + l) = ^t{t) + A^j(t). With a probabihty 
Pf, division i does not change its size — so that 
£,i{t + 1) = £,i{t) — and an altogether new division 
j is created with size £,j{t -I- 1) = A^i(t). Thus, the 
parameter p f reflects the tendency to diversify: the 
larger is , the more likely it is that new divisions 
are created. 

The dynamics are thus controlled by four parameters: 
V, Pa, and Pf, S'inin just sets the scale, so the re- 
sults of the model do not depend on its value. We assume 
that there is a broad distribution of values of Smin in the 
system because companies in different activities will have 
different constraints. 

In Fig. |2[ we compare the predictions of the model for 
the distribution of company sizes with the empirical data 
I psf . We find similar results for a wide range of parame- 
ters: V = 0.1 -0.2, Pa = 0.01-1, andp/ = 0.1-1.0. We 
define one "year" to be £ iterations of our rules applied 
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where fi represents the average growth rate. Moreover, 
the standard deviation cri(S) is consistent with a power 
law form 
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and for US manufacturing companies, /3 w 0.2 [T^]. Fig- 
ure ^ displays p{ri S), and is quite similar in form to 
empirical results |^|. Figure |3|d compares <yi{S) with 
the empirical data of Ref. |jl^: for both, Eq. (||) holds 
with 13 = 0.17±0.03. Equations (|)-(|) allow us to scale 
the growth rate distributions for different company sizes 
(Fig. |c). 

We next address the question of the structure of a given 
company. To this end, we calculate the probability den- 
sity pi{£,i\S) to find a division of size in a company of 
size S. For the model, we find that the distribution pi 
is peaked at a maximum which scales as S". Hence, we 
make the hypothesis that pi obeys the scaling relation 
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This hypothesis is confirmed by the scahng plot of 
Fig. I^a. We find a = 0.66 ± 0.05 from plotting the aver- 
age value of against 5*. The same value of a leads to 
the best scaling plot. 
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FIG. 3. (a) Probability density of one-year growth rates 
for firms with different sizes. The values of the parameters 
are the same as in Fig. ^. The distributions are tent-shaped, 
as for the empirical data consistent with an exponential 
distribution, (b) Dependence of the standard deviation of the 
growth rates on company size. Shown are the predictions of 
the model and the empirical results. The straight line with 
slope 0.17 is a least square fit to the predictions of the model, 
(c) Probability density of one-year growth rates for different 
company sizes plotted in scaled variables. 

Next, we make the hypothesis that the probability den- 
sity y02(^|'S') to find a company with size S composed of 
N divisions, obeys the scaling relation 
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In writing (||), we use the fact that from (Q) the character- 
istic size of a typical division scales as S" , so that the typ- 
ical number of divisions in a company is S/S" ^ S^^". 
Figure ^ shows that the results of the model are consis- 
tent with the scaling relation (||) , with the same value of 
the scaling exponent a used in Fig. |4|a. 



The results described by Eqs 
agreement with empirical studies 




are in qualitative 
that show larger 



firms to be more diversified, 
to the simple scaling law 
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which can be tested. For a = 0.66 ± 0.05, (g) predicts 
/3 = 0.17 ± 0.03, in remarkable agreement with our inde- 
pendent calculation of f3. 
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FIG. 4. (a) "Data collapse" of the conditional probability 
density pi; the data fall onto a single curve corresponding to 
the scaling form (^. (b) "Data collapse" of the conditional 
probability density p2\ the data fall onto a single curve corre- 
sponding to the scaling form (^. 

We find that the predictions of the model are only 
weakly sensitive to the parameter values, which perhaps 
is the reason why firms operating in quite different indus- 
tries are described by very similar empirical laws. Ac- 
cordingly, we conjecture that the scaling laws found for 
US manufacturing firms also hold for other countries, 
such as Japan, with (3 « 0.2; this conjecture is currently 
being tested with empirical data [ pTf . 

The present model rests on a small number of assump- 
tions. The two key assumptions are: (i) Firms tend 
to organize themselves into multiple divisions once they 
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achieve a certain size. This assumption holds for many 
modern corporations (ii) Growth rates of differ- 

ent divisions are independent of one another. For an 
economist, the latter is perhaps the stronger of the two 
assumptions. We find that correlations in the growth 
rates of divisions within a same company, even weak cor- 
relations, lead to /3 « 0. Thus, we confirm that it is the 
assumption of independence among the growth rates that 
reproduces the empirical findings of Refs. [ p^ . 

There are two features of our results that are perhaps 
surprising. First, although firms in our model consist of 
independent divisions, we do not find (3 ~ \/2. To un- 
derstand why /3 < 1/2, suppose that the distribution of 
Sm = InS'min IS & Dirac 5-function. Although this as- 
sumption is unrealistic, it leads to an understanding of 
the underlying mechanisms in the model. For this case, 
we find (i) that the distribution of company sizes is still 
close to log-normal, with a width W which is a function 
of the parameters of the model and, (ii) that the num- 
ber of divisions increases linearly with size, so a = and 
(3 = 1/2. Then, by integration over Sm, we can estimate 
the value of (3 for the case of a broader distribution of 
Sfn- Suppose that Sm follows some arbitrary distribution 
with width T>. Averaging crj(5) over this distribution, 
we find (3 = W/2(I? + yV). For a wide range of the values 
of the model's parameters, V > W, and we find that 13 
is remarkably close to the empirical value /3 « 0.2. 

Second, the distribution p[ri\S) is not Gaussian but 
"tent" shaped. We find this result arises from the inte- 
gration of nearly-Gaussian distributions of the growth 
rates over the distribution of ^min- For large values 
of |ri|, the saddle point approximation gives p{ri\S) ~ 
exp(— log^ \ which decays slower than exponentially, 
in qualitative agreement with the model's predictions and 
with empirical observations. For |ri| <C 1, p{ri\S) is ap- 
proximately Gaussian, while for intermediate values of 
|ri|, the distribution decays exponentially. Our analyti- 
cal predictions are in agreement with the model and with 
empirical results. 

The model leads to a number of conclusions. First, it 
suggests the deviations in the empirical data from pre- 
dictions of the random multiplicative process may be ex- 
plained (i) by the diversification of firms, i.e., firms are 
made up of interacting subunits; and (ii) by the fact that 
different industries have different underlying scales, i.e., 
there is a broad distribution of minimum scales for the 
survival of a unit (for example, a car manufacturer must 
be much larger than a software company). 

Second, the model suggests a possible explanation for 
the common occurrence of power law distributions in 
complex systems. Our results suggest that the empiri- 
cally observed power law scaling does not require some 
"critical state" of the system, but rather can arise from 
a interplay between random multiplicative growth and 
the complex structure of the units composing the sys- 



tem. Here we addressed the case in which the interac- 
tions between the units can be treated in a mean field 
way through the imposition of a minimum size for the 
subunits. We believe that more complex interactions will 
still lead to power law scaling, and that our model may 
offer a possible framework for the the study of complex 
systems. 
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